Supersymmetry in carbon nanotubes in a transverse magnetic field 
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Electron properties of carbon nanotubes in a transverse magnetic field are studied using a model 
of a massless Dirac particle on a cylinder. The problem possesses supersymmetry which protects 
low-energy states and ensures stability of the metallic behavior in arbitrarily large fields. In metallic 
tubes we find suppression of the Fermi velocity at half-filling and enhancement of the density of 
states. In semiconducting tubes the energy gap is suppressed. These features qualitatively per- 
sist (although to a smaller degree) in the presence of electron interactions. The possibilities of 
experimental observation of these effects are discussed. 



I. INTRODUCTION 

The electronic properties of single-walled carbon nan- 
otubes (NT) vary for tubes with different structure. De- 
pending on the angle between the tube axis and graphite 
lattice, called the NT chiral angle, the tube electron spec- 
trum can be metallic or semiconducting. 1 The semicon- 
ducting band gap of a single-walled NT (SWNT) is of 
the order of 1 eV and scales inversely with the NT radius 
R. Also, in many nominally metallic tubes a minigap ap- 
pears at the band center due to the intrinsic tube curva- 
ture. 2-4 These gaps, recently observed experimentally, 5,6 
have the size of a few tens of millivolts for SWNT's and 
scale as 1/R 2 with the NT radius. 

Ajiki and Ando have made a remarkable observation 7 
that in nanotubes the metallic behavior is fragile: a 
metallic NT can be easily turned into a semiconducting 
one by applying a relatively weak parallel magnetic field. 
Such a field, by inducing backscattering between right 
and left electron modes, opens a minigap at the band 
center. This gap, linear in the field, is given by the mag- 
netic flux scaled by the flux quantum, ttR 2 B/& , times 
the semiconducting gap size. For B ~ 10 T the gap is 
of the order of lOmeV for typical SWNT radii. Effects 
of parallel field on multi-walled NT's have been reported 
in Ref. 8 . Electronic properties are also sensitive to me- 
chanical distortion, such as twisting, bending, or squash- 
ing, 9-11 as well as to external electric fields. 12,13 

Another interesting observation made in Ref. 7 is that a 
transverse magnetic field affects electron states in a way 
completely opposite to the parallel field effect. In metal- 
lic tubes the Fermi velocity is suppressed by a transverse 
field, while the density of states near the band center 
is enhanced. At the same time, in semiconducting tubes 
the band gap is suppressed. The goal of the present work 
is to rationalize these properties using the concept of su- 
persymmetry. 14 Supersymmetry has a profound effect on 
the low-energy properties by protecting the states at the 
band center. We derive a supersymmetric Hamiltonian 
for and present a simple analytic theory of the above ef- 



fects. We find that the metallic behavior is protected 
by the supersymmetry for any magnetic field that is ap- 
plied perpendicularly to the NT and does not vary along 
the tube. For a uniform transverse field, in particular, 
the Fermi velocity in metallic tubes is suppressed by the 
uniform transverse field as 

v(B)=v/I (x), x = 4nR 2 B/<I> 0l (1) 

where v is the Fermi velocity for B = and Iq(x) is the 
modified Bessel function. The density of states near the 
band center is enhanced by the same factor Iq(x). We 
also calculate the suppression of the gap in semiconduct- 
ing tubes. 

The typical field strength required to make these effects 
pronounced is quite high. The fields necessary to sig- 
nificantly alter the electron dispersion can be estimated 
from irR 2 B ~ $o = hc/e, which gives many tens of tesla 
even for the tubes of the largest available radii. These 
fields are not hopelessly strong — they are available, for 
example, in pulsed magnetic field sources 15 which allow 
one to reach fields up to 100 T. 16 Because of that and 
also because of the novel features arising from the super- 
symmetry and Dirac character of low-energy states, we 
believe that this problem is sufficiently interesting. 

The paper is organized as follows. We first review the 
basics of the carbon 7r-clcctron tight-binding band and 
its relation with the massless Dirac equation, paying par- 
ticular attention to coupling to external fields. Then we 
present a theory of a massless Dirac particle on a cylinder 
in a transverse magnetic field and calculate the spectrum 
and density of states. We analyze both the metallic and 
semiconducting tubes. Then we briefly discuss the be- 
havior in extremely high fields, where a connection can 
be drawn with the Landau levels and snake states 17 con- 
sidered previously in the context of the quantum Hall 
effect. 

After that we discuss the effects beyond the Dirac ap- 
proximation arising from the next order in the gradient 
expansion of the tight-binding problem. These effects are 
small but interesting, because they violate supersymme- 
try and lead to minigaps appearing at the band center. 
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The effects beyond the Dirac model are controlled by the 
so-called trigonal warping interaction. We consider it in 
the presence of a magnetic field and show that its effect 
depends on the chiral angle and, in particular, is absent 
for truly metallic armchair nanotubes. These effects have 
been discussed in Refs. 18 ' 19 , using a combination of nu- 
merical and analytic methods, for zigzag and armchair 
tubes. We extend the results of Refs. 18 ' 19 by consider- 
ing nanotubes with arbitrary chiral angle and also in the 
presence of minigaps of other origin. 

We also discuss the experimental implications of this 
work. The gap suppression in semiconducting NT's wins 
over the Zeeman splitting at reasonable fields. We con- 
sider the possibility to observe the gap suppression in the 
tunneling density of states and in the thermally activated 
transport regime. 

Finally, we consider the competition between super- 
symmetry and electron-electron interactions. We find 
that although the effect of supersymmetry is reduced by 
strong interactions, the qualitative features of the spec- 
trum (suppression of the plasmon velocity and of the 
scmimetallic gap) are similar to those of the noninter- 
acting case. 



II. DIRAC MODEL FOR THE CARBON n 
BAND AND NANOTUBES IN EXTERNAL 
FIELDS 

Here we review the basics of the theory of electron 
states of the two-dimensional (2D) carbon monolayer, 
making a connection with the 2D Dirac equation. This 
will provide a good starting point for the following dis- 
cussion of nanotubes in external fields. We shall start 
with the tight-binding description of the carbon ir band, 
following the approach of Ref. 20 , recall how the Dirac 
equation arises in this system, and then consider elec- 
tron coupling to external electromagnetic fields. 

The tight-binding Hamiltonian on a honeycomb lat- 
tice of carbon atoms with hopping amplitude t between 
adjacent sites has the form 



&l>{r) = -t J2 ^ r ') 



(2) 



\r'— r\=a c 



where r' are the nearest neighbors of the site r, and 
a cc is interatomic spacing. In carbon, t w 3 eV and 
a cc — 0.1437 nm. For simplicity and because the elec- 
tron spectrum is e — > — e symmetric, from now on we 
shall ignore the minus sign in Eq. (2). 

The zero chemical potential in Eq. (2) describes the 
half- filled tt band, i.e., the density of one electron per 
site. For an infinite system, the states of the prob- 
lem (2) are plane waves and the spectrum is given by 
e(k) = ±t| e lk ■ r *|, where are the nearest- neighbor 
bond vectors. This is a spectrum of a semimetal with 
the conduction [e(k) > 0] and valence [e(k) < 0] sub- 
bands touching each other at two points K and K' in 
the Brillouin zone. 



The tight-binding bandwidth 6i ~ 18 eV is much larger 
than the energies of the states close to the band center 
considered below. Because of that, it is useful to project 
the problem (2) onto the subspace of states with |e| <C t 
and derive an effective low energy Hamiltonian for such 
states. To carry out the projection, we note that there 
are only four independent states with e = 0. These states 
form two complex valued conjugate pairs which we de- 
note as u(r), v(r) and u(r), v(r). It is convenient to 
choose the states u and v to be zero on one of the two 
sublattices of the honeycomb lattice. On the other sub- 
lattice each state takes the values 1, u> = e^ 27r / 3 ^ and 
CD = e~( 27r / 3 ) J (see Fig. 1). The states u{r) and v(r) have 
the same quasimomentum of a value 



K 



A-k 



(3) 



opposite to that of the states u(r) and v(r). Each pair 
of states u(r), v(r) and u(r), v(r) forms a basis at the 
points K and K 1 , respectively. 

Projecting the wave function ip(r) on u(r) and v(r) 
and, respectively, on u(r) and v(r) defines Dirac spinor 
components for each of the two points K and K' . We 
focus on the u, v pair and write the states near the point 
K with small energies |e| <C t as linear combinations 



ip(r) = ipi(r)u(r) + ip 2 (r)v(r), 



(4) 



with the envelope functions ipi,2{ r ) varying on the scale 
much larger than the interatomic spacing a cc . By substi- 
tuting the wave function (4) in the tight-binding Hamil- 
tonian (2) we have 

eipi(r) = t{ip 2 (r—a) +w^(r-wa) + w^(r-wa)} , (5) 
ei/) 2 (r) = t{ipi(r+a) + Luipi(r+Lua) + w^i(r+wo)} , (6) 

where a is a shorthand notation for a cc . Here the prod- 
ucts za with unimodular complex numbers z = 1, u), CD 
in the arguments of ipi,2 are understood in terms of 2D 
rotations of the vector ax by argz. 

Expanding slowly varying ipi.2(r), we obtain 



eV>i(r) = -hv(d x - id v )^ 2 {r), 
etp 2 (r) = hv(d x + idy)ipi{r), 



(7) 



where v 



\t a cc /h. The Hamiltonian (7) defines mass- 



less Dirac fermions with the linear spectrum e(fc) = 
±?w|k|. In carbon, the velocity v = 8 x 10 7 cm/s. Similar 
relations hold for the point K' . 

Equations (7) can be cast in the conventional Dirac 
form eip = Hip with 



H = v a ■ p = v (aipi + a 2 p 2 ) 



(8) 



for the two-component wave function ip = (ipi,ip 2 ) T , with 
Q!i,2 given by the Pauli matrices: 



ai = a 2 , a 2 



(9) 
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FIG. 1: Shown are two plane-wave basis states u(r) and v(r) 
of the problem (2) with e = 0. Both u(r) and v(r) take the 
values 1, lj = e (2?r/3)i and uj = e " (2,r/3)i on one sublattice 
and vanish on the other sublattice of the honeycomb lattice. 
The states u(r) and v(r) have the same quasimomentum and 
form a basis of the Dirac problem (9) at the point K. The 
independent basis states at the point K' are u(r) and v(r). 

The Hamiltonian near the K' point can be derived in 
a similar way. The result has the form (8) with a sign 
change in the second term: u\ — a 2 , ce 2 = o\. 

Below we shall consider electrons in the presence of 
external electromagnetic fields. The minimal form of the 
coupling to external fields follows from the gauge invari- 
ance: 

H = va ■ (p - -A) + eip , (10) 

where ip and A are the scalar and vector electromagnetic 
potentials. The effect of electron spin, ignored here for 
simplicity, can be included in Eq. (10) via a Zeeman term. 

Equation (10) describes the lowest-order approxima- 
tion in the gradients of ^1,2 and the potentials ip and A. 
Here we consider the exact tight-binding equations in the 
presence of external fields: 

eVi(r)=i( Y, ze l ^—^ 2 {r-za)\, (11) 

\Z — 1,UJ,UJ / 

eip 2 ( r )=ti ze^r+^iir + za) j , (12) 

\Z — 1,CJ,LJ / 

where the phases j r)r / are the integrals of the vector po- 
tential along the nearest-neighbor bonds, 

r 

7 = Jl f A (x) • dl . (13) 

r' 

Equations (11) and (12) can be used to obtain the gra- 
dient terms of higher order along with the coupling to 
external fields. One can check that expanding the ex- 
ponents in Eqs. (11) and (12) and keeping the lowest 
nonvanishing terms gives the Dirac Hamiltonian (10). In 
Sec. VI we shall use Eqs. (11) and (12) to obtain higher- 
order corrections to Eq. (10). 



To apply the above results to nanotubes, we consider 
electrons on a carbon sheet rolled into a cylinder. The 
transformation of the tight-binding problem (2) to the 
Dirac problem (7) based on the representation (4) is valid 
provided that the cylinder circumference L — 2ttR is 
much larger than the interatomic spacing a cc . Since for 
typical NT radii the ratio L/a cc can be between 10 and 
20, the approximation (4) is entirely adequate. 

The NT electron properties, depending on the nan- 
otube structure, can be either metal like or dielectric 
like. Which of these situations takes place depends on the 
manner the cylinder is obtained from the carbon mono- 
layer. In the Dirac approach, the condition for metallic 
behavior can be formulated directly in terms of the func- 
tions u(r) and v(r): The nanotube is metallic if and only 
if one can define on the NT cylinder the two functions 
u(r) and v(r) according to Fig. 1 without running into 
a mismatch of the function values upon the cylinder clo- 
sure. 

To demonstrate this, let us suppose that the func- 
tions u(r) and v(r) on the cylinder exist. Without loss 
of generality we choose the x axis along the cylinder 
and the y axis along the circumference. The problem 
(7) has periodic boundary conditions in the y direc- 
tion, and thus the wave functions can be factorized as 
i/j 12 (r) = Tpi,2(x)e lknV , where k n = 2mjL = n/R. Then 
the dispersion relation for the ID problems describing 
motion along the x axis with fixed k n is 

e n {k x ) = ±hv{kl + kl) 1 ' 2 . (14) 

In this case the subband with n = has metallic proper- 
ties and the subbands with are dielectric. 

Now let us consider the other possibility when the 
cylinder is constructed in such a way that the func- 
tions u(r) and v(r) cannot be defined without a value 
mismatch. In this case, upon rolling the carbon sheet 
into a cylinder, the sites with different function values 
shown in Fig. 1 are glued together. However, since 
all values of the functions u(r) and v(r) are powers of 
lu = e( 2,r / 3 ) 1 , one notes that Eqs. (7) can still be used here 
if they are augmented with quasiperiodic boundary con- 
ditions, 4>i, 2 (x,y + L) = uip lt2 (x,y) or ip 1<2 (x,y + L) = 
&' t Pi,2(%,y)j which, combined with the value mismatch 
of u(r) and v(r), make ip(r) single valued. Factoring 
the wave function as above, one obtains ID subbands 
with the dispersion of the form (14), in this case with 
k n = (n ± Note that in this case all spectral 

branches have dielectric character. 

Now we consider a nanotube in the presence of a paral- 
lel external magnetic field. In this case, electron proper- 
ties are described by the Dirac equation (10) with ip = 
and the vector potential A with just the y component, 
A y = where $ = ttR 2 B is the magnetic flux. The 

boundary conditions in the y direction are periodic for 
the metallic case and quasiperiodic for the dielectric case. 
In the presence of a parallel magnetic field the problem 
remains separable and thus the wave function can be fac- 
torized in just the same way as above. One again finds 
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ID subbands with the spectrum (14), where 



k n R 



n± 3 + <t>\\ 



metallic, 
semiconducting, 



(15) 



for the metallic and semiconducting NT's, respectively, 
with </>|| = $/$o and $0 = hc/e. Thus in the presence of 
a parallel field the gapless n = branch of the metallic 
nanotube spectrum (14) acquires a gap. 7 ' 19 Interestingly, 
there is no threshold for this effect, since the gap forms at 
arbitrarily weak field. The gap size is 2A = 2\<fr\\ \hv/R. 
One notes that the field-induced gap appears not at the 
Fermi level but at the center of the electron band. Thus 
it affects the metallic NT properties only for electron 
density sufficiently close to half-filling. 



III. DIRAC EQUATION AND 
SUPERSYMMETRY 

At the energies smaller than the total bandwidth 6t 
(ca. 18 eV in carbon) electron states are described (sepa- 
rately near each of the K and K' points) by the massless 
Dirac equation (10). For a uniform transverse magnetic 
field, the field component normal to the NT surface is 
B±(0) = Bsin9, where 9 ~ y/R is the azimuthal angle. 
The corresponding vector potential can be chosen along 
the tube axis x, A(r) = xA(9), where 



2</> cos 9 , <j) = 



ttR 2 B 



(16) 



In this case the longitudinal momentum fik is conserved 
and the states on the NT cylinder have a plane- wave form 
ip(r) = ip(x,9) = i}j{9)e lkx . The Dirac Hamiltonian for 
i>{6) is 

U D = Aoiia^e + (kR- 2<j) cos<9)cr 2 }, (17) 



with <7i ; 2 the Pauli matrices and 

hv 

A °=R- 



(18) 



The equations near the K' point have the form (17) with 
a sign change in the first term, o\ — > —<j\. 

The eigenvalues of the operator (17) give the electron 
dispersion relation e{k). We have chosen the dimen- 
sionless transverse field parameter (f> in the form (16), 
which makes contact with the parallel field problem, 7 ' 19 
Eq. (15). 

The NT states are described by quasiperiodic wave 
functions on the cylinder (y = R9), 



i>(y + L) 



2ni8 



2irR 



with 



5 = 



, metallic, 

±| semiconducting, 



(19) 



(20) 



We consider the problem (17) with an arbitrary phase d 
in the boundary conditions (19). This will permit us to 
generalize the results to the cases of metallic NT's with 
a minigap induced by curvature 21 or in the presence of a 
parallel magnetic field. These problems can be described 
using the boundary conditions (19) with 5 slightly shifted 
away from the ideal values (20). 

The electron bands e(k) can be studied using the trans- 
fer matrix. We integrate the Dirac equation in the inter- 
val < 9 < 2ir and write a formal solution \p{9) to the 
problem H D tp = ^ as ip(9) = S(9)ip(0) with the 2 x 2 
matrix 

S(9) = Texp / {-ieat + (20 cos 9' - kR)a 3 }d9' , 

(21) 

where e = e/Ao. Here T stands for operator ordering 
with respect to 9. The quasiperiodic boundary condition 
(19) requires 



tr 5, 



0=2tt 



2 cos(2tt5) 



(22) 



Different energy bands e = e„(fc) can be found numeri- 
cally as solutions of Eq. (22). The bands obtained in this 
way are displayed in Fig. 2. 

In the present section we show that supersymmctry 
allows one to make rather general statements about 
the low-energy NT spectrum. Originally, supersymme- 
try was suggested 22 as a special symmetry between the 
bosonic and fermionic sectors of relativistic field theo- 
ries that protects the zero-energy eigenstate. Later, the 
concept of supersymmetry was brought to single-particle 
quantum mechanics. 23 This has yielded the classification 
of exactly solvable potentials using factorization of the 
Schrodinger equation (see Ref. 14 for a review). 

Below we apply the arguments of supersymmetry to a 
problem of NT electrons in a generic magnetic field B(r) 
that (i) is perpendicular to the NT axis and (ii) does 
not vary along the NT. In this case the field component 
normal to the NT surface is a function only of y = R9, 
B± = B(y). Similarly to Eq. (16), we choose the follow- 
ing gauge: 

A. = AW = %%, A y = 0, (23) 



where 



The function ip(y) is uniquely defined by demanding pe- 
riodicity 



<P(V + L)= ip(y) 



and zero average 



/ 

Jo 



dyip = 



(25) 



(26) 
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Wave vector kR 



FIG. 2: Electron dispersion e(fc) in the presence of a large 
uniform transverse magnetic field (f> = 0.5 (bold lines) and 
in the absence of the field (dashed lines). Top: metallic NT 
(5 = 0). No gap opens due to supersymmetry, but the velocity 
is suppressed according to Eq. (45). Bottom: semiconducting 
case (8 = 1/3). The energy gap is suppressed by the factor 
!?i/3(0); see Eq. (58) and Fig. 4. 



The corresponding Dirac Hamiltonian reads 



Hi) : : //," { />] (),, ' { k - ^ J 17 -2 



(27) 



which reduces to Eq. (17) when B(y) = Bsm(y/R). It is 
useful to decompose Hd into two pieces, 



Hd = 



)t 



where 



and 



Q 




A 



A — hv { id v + i I k — 

dy 



A* 




dip 



(28) 
(29) 



hv < id. 



dtp 
dy 



(30) 



The connection with the supersymmetric quantum me- 
chanics (see Chap. 2 in Ref. 14 ) is established by con- 
structing a supersymmetric Hamiltonian 



HsUSY 



(J$A 
A4* 



(31) 



which, together with Q and Q\ satisfies the superalgebra 
41/1), 

[HsUSY, Q] = [HsUSY, Q t ] = 0, /™i 

{Q,Qt} = Hsusy, {Q 1 Q} = {Q\Q^} = 0. 

Here {A, B} = AB + BA stands for anticommutator of 
operators. In relativistic field theories, the supercharges 
Q and transform fermionic and bosonic degrees of 
freedom into each other. Although we deal with a sin- 
gle electron NT spectrum, one could formally interpret 
the upper and lower components of the wave function ip 
as fermionic and bosonic sectors of the supersymmetric 
Hamiltonian Hsusy- 

One interesting implication of the algebra (32) is that 
Hsusy can be expressed as sums of the square of Hermi- 
tian supercharges, Q\ and Q 2 , 



USY 



— Qi + Qi, 



where 



Q 1 = _(Q + Qt) 



Q-2 



From Eq. (29), one verifies Q\ = Q\ and 
Hsusy = 2Q 1 = H 2 D - 



(33) 



(34) 



(35) 



Thus the energy spectra of Hsusy and Hd are closely 
related. 

Let us now show how supersymmetry protects zero- 
energy states of Hsusy and Hd- For that, we construct 
zero energy states of Hsusy that are compatible with the 
boundary condition (19). Due to Eq. (33), any such state 
tp satisfies 



Qiip = Q2V' = 

or, equivalently, 

Qtp = = . 
The latter equation has two independent solutions 



— p-ky+v(y) 



lj'2 



(36) 



(37) 



(38) 



Note that since k is real and <p(y) is periodic in y, 
the zero-energy solutions (38) are compatible with the 
boundary condition (19) if and only if S = and k = 0. 
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For the latter case the exact zero-energy eigenstates of 
Tijj can be written as 

where normalization requires 

LgP = [ dye 2 ^ , Lg (2) = C dye~ 2 ^ . (40) 
Jo Jo 

In the case of a uniform perpendicular field 

p = 20 sin(9 , 9 = y/R, (41) 

the normalization factors are given by the modified Bessel 
function: 

.9o 1} = .9c- 2) = ^ / e**«* 9 dB = h<M) . (42) 

The states (39) are degenerate at any field strength. 

Let us stress here that the zero-energy eigenstates (39) 
of Hd exist 24 at 6 = 0, k = for a generic magnetic 
field perpendicular to the NT and not varying along the 
tube. An example is a field of a current flowing along 
the wire parallel to the nanotube axis. Moreover, the 
above arguments still apply if the NT does not have a 
circular cross section, as long as the minigap due to the 
curvature 2,3 of the graphene sheet is not open. In what 
follows we confine ourselves to the case of a cylindrical 
NT in a homogeneous perpendicular magnetic field for 
simplicity, bearing in mind the generalizations mentioned 
above. 

Using the states (39) one can study how the linear 
dispersion relation changes near the band center e = 0. 
For that we project the Hamiltonian (17) onto the basis 

i> 1<2 (x,e) = e ik *^l(9) . (43) 

The projected Hamiltonian 

_ v 
Hd * m = hkv (T 2 , v = (44) 

yields the dispersion relation 

e(fc) = ±M (ioM) • (45) 

This describes a reduction of the Fermi velocity h~ 1 de/dk 
near e = by a factor 7 (4</>). Since Io(i<j)) > 1, the 
density of states at the band center, 

v = dN/de=^I (4<f>), (46) 

is enhanced. (The factor of 4 accounts for the spin and 
valley degeneracy neglecting the Zeeman splitting; see 
Sec. IV below.) Due to the exponential behavior of the 
Bessel function in Eq. (45) at large c/>, this enhancement 




-2-1012 

Energy e/A 



FIG. 3: Density of states v — dN/de in metallic NT's per one 
fermion flavor as a function of energy in the units of Ao = 
hv/R. The peak value v(0) at the band center is given by 
Eq. (46). 

becomes dramatic at high fields (Fig. 3). The reduction 
of the Fermi velocity and the corresponding density-of- 
states enhancement is a general consequence of a super- 
symmetry. Indeed, in the generic field case described 

above the velocity is reduced by a factor *J g\p gip > 1. 

Both g^\ g^ 1 > 1 [see Eq. (40)] since the average of an 
exponential is greater than or equal to an exponential of 
the average. 

In the single-particle approximation, the tunneling 
density of states coincides with the thermodynamic den- 
sity of states (46). The peak in the tunneling density 
of states at the band center is a distinct manifestation 
of supersymmetry. Recently, a scanning tunneling probe 
has been used 6 to study curvature-induced minigaps in 
nominally metallic tubes placed on a metallic substrate. 
In this setup, the electron interactions that could have 
modified the single-particle behavior are screened by the 
substrate, and the measured density of states is unaf- 
fected by Luttinger liquid effects. In a similar system in 
a high transverse field, an enhancement of tunneling at 
the band center, Eq. (46), and a peak in the density of 
states (Fig. 3) are the experimental signatures to look 
for. 



IV. GAP SUPPRESSION 

In the present section we consider nanotubes that have 
a gap A at the band center. We characterize them by the 
parameter \8\ — A/A that enters the boundary condi- 
tion (19). There are two kinds of such tubes: semicon- 
ducting NT's with 5 = ±1/3 and nominally metallic NT's 
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in which a minigap appears due to curvature or external 
field, yielding a small \S\ <C 1. The Dirac problem (17) 
and (19) is not supersymmetric for 5^0. However, since 
the supersymmetry is an exact property at S = 0, one 
can still expect it to manifest itself in a problem with 
a relatively small S. Below we show that a gap at the 
band center is suppressed in the presence of an external 
transverse field: 

A W = ^T' 9(W>1, (47) 

gs{4>) 

where the gap suppression factor g$((j>) diverges as <j> — > 
oo. This means that supersymmetry is restored in the 
limit of a strong field. The suppression of the gap in 
semiconducting NT's is illustrated in Fig. 2 (lower panel). 

We shall first consider a simpler case of a nominally 
metallic tube with j « 1. The gap in this case can 
be found using perturbation theory in 5. For that, we 



perform a gauge transformation ij)(ff) 



ip'(6) and 



for if)'(d) obtain a problem with the periodic boundary 
condition. The new Hamiltonian differs from Eq. (17) by 
a term linear in 5: 



H' = H D - £A cn . 
It is convenient to rewrite this Hamiltonian as 
H' = H$ + A (fci?CT2 — 5<j\) , 



(48) 
(49) 

where Tt'ff is the Hamiltonian (17) with k = 0. We 

note that is a supersymmetric Hamiltonian with 

eigenstates (39). The spectrum of the Hamiltonian Ti' 
at small kR and 5 can be found by projecting the second 
term of Eq. (49) on the basis (43) of plane-wave states 
constructed out of Eq. (39). This yields the dispersion 
relation 

Thus we find that in this case the gap is suppressed by 
the same factor (42), 



<7o(0) = /o(40) > 1, 



(51) 



as the Fermi velocity in metallic NT's. The gap suppres- 
sion is described by Eq. (51) in the limit of small 8 for 
any magnetic field 0. 

One can also study the gap suppression analytically 
for generic S using perturbation theory in the field 0. 
The energy e(k = 0) that gives the gap is defined by the 
condition (22) for the transfer matrix. We calculate the 
trace of the transfer matrix (21) at k = by expanding 
it perturbatively in </> ^ 1: 

tr S e =2* = 2cos27re + 8A 1 2 + 32A 2 4 + e>(0 6 ), (52) 

where e = e(0)/A o and the coefficients Ai i2 are given by 

27re sin 27T6 

Ai = ; — , (53) 



A 2 = -7re 



1 - 4e 2 ' 

2?re (l - 4e 2 )cos27re+ (§ + 6e 2 )sin 2vre 
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FIG. 4: Gap suppression factors (47) for nominally metallic 
NT's, 5 = (solid line), and semiconducting NT's, 5=1/3 
(dashed line), as a function of field <f>. The fine solid line is 
the result of the expansion (58) in powers of <j>. (Note the 
logarithmic scale for gs.) 



The condition (22) on the energy along with the defini- 
tion of the suppression factor (47) gives 

gs = 1 + a 5 (t, 2 + + O^ 6 ) , 

4 a 1 - 20£ 2 

Ps 



(55) 
(56) 



1-4<S 2 ' {1-4S 2 Y 
Substituting, in Eq. (56), 6 — and 5=1/3 we obtain 
go ((/)) = l + 4(f> 2 + 4(/> 4 + 0(c/) 6 ) (metallic) , (57) 

5i/3(0) 



1 

3fi 396 
1 + ^</> 2 - + O(0 6 ) (semiconducti( 

5 25 



~2\3 



;i-4?) 



(54) 



The expression (57) coincides with the Taylor expansion 
of I (40). 

These analytical results can be compared with the gap 
suppression factors obtained numerically (Fig. 4). For 
nominally metallic NT's with small minigap we find that 
at S <C 1 the value gs~*o is accurately given by Eq. (51). 
The analytical expression (51) coincides with the numer- 
ics in the entire field range. In the semiconducting case 
of S = 1/3 the expansion (58) works reasonably well at 
< 1/4. At larger fields (f> > 1/4 the gap is suppressed 
exponentially, 51/3(0) oc e 4 ^ (see Fig. 4). 

Let us discuss the possibilities to observe the suppres- 
sion of the gap. A competing effect due to the magnetic 
field that leads to a gap suppression is the Zeeman spin 
coupling Hz = —jJLoB with fi = eh/2mc. The gap sup- 
pression at weak fields, <C 1, with the Zeeman effect 
included, is described by 

A(0) = A - a«50 2 A - /jtB, (59) 
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with as given by Eq.(56). The Zeeman effect, linear in 
B, dominates at weak fields. However, the orbital effect 
as<j) 2 A, quadratic in B, overcomes the Zeeman effect at 
relatively moderate fields <f> <C 1. 

For semiconducting NT's the gap A = hv/3R, and 
Eq. (56) gives ai/ 3 = 36/5. In this case, the inequality 
cui/3(/> 2 A > \xB yields 4> > 1.13 cib/R with as = h 2 /me 2 
the Bohr radius. Using carbon parameters we estimate 
that the magnetic field has to exceed 

B = 78(i?[nm])- 3 T (60) 

which can be low enough for tubes of large radius. 

The energy gap can be studied experimentally by mea- 
suring transport in the thermally activated regime. The 
thermal activation energy will depend on the magnetic 
field according to Eq. (59). The resistance 

R(T) cx exp[A(0)/fc B T] (61) 

will be sensitive to magnetic field because the variation of 
the gap can exceed k B T even at the fields much smaller 
than Eq. (60). For example, consider a NT of radius 
R = lnm, in which case A = hv/3R = 0.178 eV. For 
the magnetic field B = B from Eq. (60) we have <f) = 
0.060 and ai/ 3 (f> 2 — 0.026. In this case the gap shift 
a 1/3 (/) 2 A = 4.6 meV is larger than k B T at T < 53 K. 

V. ELECTRON SPECTRUM IN EXTREMELY 
LARGE FIELDS 

Below we consider the qualitative features of the en- 
ergy bands in metallic NT's in the limit of a large uniform 
external field, corresponding to <fi > 1. The behavior of 
the band dispersion at \k\R <C 4> can be attributed to 
the Landau levels of the problem (17). Let us consider a 
square of the Hamiltonian (17) and the eigenvalue prob- 
lem Ttjjip — e 2 ip. It reads 

H 2 D i> = A 2 [-d 2 + U k {9)]^ = e 2 W V> , (62) 
U k (9) = {2cj) cos 6-kRf - 2<jxr 3 sin 9 , (63) 

with ip(9) obeying the boundary conditions (19). Note 
that Hjj is a diagonal 2x2 matrix in the space of spinors 
tp. In what follows we take 03 = +1. 

When cj) ^> 1, the kinetic energy of Eq. (62) is frozen 
since Uk{0) cx 2 is much greater than dg. In this limit, 
the Hamiltonian is dominated by the potential energy 
term U k (9) and the low-energy states are localized near 
the minima of U k (9). At <j> 3> 1 the potential U k (9) 
has two slightly asymmetric minima near 9± — ±7r/2, 
where it can be approximated by a harmonic potential 
U k (9) w =f20 + 4</> 2 (# - 0±) 2 - The size of electronic wave 
function in the circumferential direction is 

w~l B =R/^/24><^R, (64) 

where Ib = (hc/eB) 1 / 2 is the magnetic length. Thus at 
large field the electrons are localized near the extrema 



of the magnetic field 9±. In these regions the field is 
approximately constant: |_Bj_(#)| w B. 

The Landau level spectrum e 2 (fc) of Eq. (62) obtained 
within the harmonic approximation yields approximately 
^-independent levels for Hd- 

e n {kR^4>) = ±2A 0V / n\^\ , n = 0,1,2,... . (65) 

The lowest-energy level of Eq. (62) is e = at k = 0. This 
value, because of the supersymmetry, is not approximate 
but exact. 

The behavior at large momenta |fc|i? 4> can be un- 
derstood semiclassically in terms of the so-called snake 
states. 17 Snake states correspond to a classical particle 
moving along zero-field lines. This motion is stable for 
a particle traveling in one direction and is unstable for 
it traveling in the opposite one. The snake states are lo- 
cated at 9 = and 9 — ir, where the field B{9) = Bsin9 
vanishes. This is consistent with the high-field limit of 
the problem (62), since at \k\R 3> 4> the minimum of 
U k (9) is 9 w for kB > and 9 w tt for kB < 0. The 
dispersion relation for such states is 

e(fc) - ±Aoy/U k (0o) w ±hv(\k\ - 2\4>\/R) . (66) 

This linear dispersion relation with an offset 2\(f>\ Uv/R 
holds even for small fields as long as \k\R ^> 
max {<fi, l/(j>} (see Fig. 2, large kR). 

Since # is different for positive and negative k, the 
left- and right-moving snake states are spatially sepa- 
rated. For 4> > 0, for instance, the left- (right-) moving 
snake states are localized near 9 = tt (9 = 0). The char- 
acteristic width of the snake-state wave function is 

^snakc = ^ /2 (i?/fc) 1/4 OC ^B^ 1 ^ . (67) 

The width w sna k c <C R when kR ^> 1/|0|. In the Dirac 
problem the wave function width (67) is different from 
that for the Schrodinger problem discussed in Ref. 17 , 
where w snak o oc |VB|~ 1/2 . 

VI. BEYOND THE DIRAC EQUATION: 
SUPERSYMMETRY BREAKING 

The two effects considered in Sees. Ill and IV — the 
density of states enhancement at e = in metallic NT's 
and the suppression of the energy gap in semiconducting 
NT's — are manifestations of the supersymmetry of the 
low-energy Dirac Hamiltonian (17). However, supersym- 
metry is not present in the original tight-binding problem 
(2). Below we show that the terms correcting Eq. (17) in 
the next order of the gradient expansion violate super- 
symmetry. Thus the supersymmetry in nanotubes is not 
exact but approximate: the nonsupersymmetric effects 
are small in a cc /R. 

To obtain the nonsupersymmetric terms of the effec- 
tive Hamiltonian we consider the low-energy subspace of 
states with |e| <§C t near the Dirac points K and K' . The 
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basis states at e = are the functions u(r), v(r) and 
u(r), v(r) defined in Sec. II (see Fig. I). The wave func- 
tion near the point K (K 1 ) can be represented as a lin- 
ear superposition (4) of u(r) and v(r) [respectively, u(r) 
and v(r)] multiplied by the smooth envelope functions 
ip 1,2 (r). We perform a gradient expansion of the slowly 
varying envelope functions using Eq. (11). In the lowest 
nonvanishing order in a cc dipi,2 we retain the Hamiltonian 
(17) with v = 2>ta cc /2h. 

The terms of second order in the gradients a1 c didjipi,2 
give the required correction H tvJ to the Dirac Hamilto- 
nian Hd called the carbon trigonal warping interaction. 
In this case, since we are interested in the problem in an 
external field, the gradient expansion of Eq. (11) should 
be accompanied by an expansion of the phase factors 
(13). After this expansion is carried out we choose the 
tube axis orientation with respect to the carbon lattice 
by specifying the chiral angle 6. The full Hamiltonian 
Htot obtained in such a way for NT's has the form 



>(e/2)<r 3 



-10CT3 



H tw e 



iO<T3 



(68) 



where Hd is the Dirac Hamiltonian (17) and the trigonal 
warping interaction Htw is given by 

Wtw = -|g Ao {(« 2 + dl) ai+i (iKde + ^j 02 



where 



k(9) = kR-^A(0) 

he 



(69) 



(70) 



The term Htw breaks the supersymmetry of the Hamil- 
tonian Wtot- Thus we expect the zero-energy state to 
disappear. Note that 7i tw also breaks the rotational sym- 
metry of Htot since 9 cannot be removed from Htot via a 
unitary transformation. Thus the behavior of the energy 
gap in a transverse field will in general depend on 6. It 
can be verified that, in the absence of external fields, the 
energy spectrum of Htot is periodic in 9 with period 7r/3, 
which is a manifestation of the 60° rotation symmetry in 
the honeycomb lattice. 

It is explicit in Eq. (69) that the effects of Htw are the 
order of a cc /R corrections to Hd- Such effects are neg- 
ligible for semiconducting NT's because of a large gap 
Ao/3. In metallic NT's, however, the Hamiltonian Htw 
plays an important role. In particular, the system de- 
scribed by the Hamiltonian (68) develops a minigap 



2a cc 4> 2 1 cos 39 1 



9o{4>) 



(71) 



due to the magnetic field <j> [defined in Eq. (16)]. This 
result can be obtained by projecting the perturbation 
Htw taken at k = onto the supersymmetric basis (39): 



H 



2a c 



tw|fc=0 — ~~ " 



R 



■An 



5o (0) 



(72) 



Note that the minigap (71) depends explicitly on the chi- 
ral angle 9 as | cos 39 1. Thus, for a given NT radius, A tw 
reaches its maximum in zigzag NT's (9 = 0) and van- 
ishes in armchair NT's (9 = 7r/2). The gap (71) is a 
manifestation of the broken supersymmetry. Minigaps of 
purely magnetic origin have been reported in Ref. 19 for 
zigzag and armchair NT's; however, the gap dependence 
on the chiral angle 9 was not discussed. 

When the magnetic field is large, <f> ~ 1, the minigap 
Atw is comparable to the curvature-induced minigap 



A c = — — A | cos39| 
16i? 1 1 



(73) 



which is present in zero field. 2,3 In this situation, the 
two mechanisms for minigaps may compete and should 
be considered simultaneously Instead of imposing the 
quasiperiodic boundary condition (19), we take the cur- 
vature effect into account in an alternative way 3 by in- 
troducing a pseudovector potential A^ c \ 



A (C) , A A (C) _ A £0 . 



„3i6 



2tt 16i? 2 



(74) 



which should be added to the magnetic vector potential 
A(r). Surprisingly, the two gap opening mechanisms, 
Eqs. (69) and (74), interfere destructively at both the K 
and K' points and produce the gap 



(75) 



ga{4>) 



in a moderate transverse magnetic field (<j> < 1). In par- 
ticular, A^, vanishes at <fi = (A^/2)" 1 due to the destruc- 
tive interference. 



VII. ELECTRON INTERACTIONS VS. SUSY 

Electron interaction effects on NT's were addressed in 
various theoretical 25-29 and experimental 30 studies. Be- 
low we consider the effects of the repulsive interaction 
between electrons on the NT spectrum in the presence of 
a perpendicular magnetic field. As we have seen above 
in Sees. Ill and IV, supersymmetry of the single-particle 
problem makes nanotubes "more metallic" in a strong 
field by enhancing the density of states and suppressing 
the excitation gap at zero energy in metallic and semi- 
conducting tubes correspondingly. It is also known that 
the repulsive interaction between NT electrons opens a 
small gap in an otherwise metallic tube 25 as well as en- 
hances the excitation gap in a semimetallic tube, 26 mak- 
ing the tubes "less metallic." Below we study the com- 
petition between supersymmetry and repulsive electron 
interactions, and find that strong interactions drastically 
reduce the effect of supersymmetry. This happens be- 
cause supersymmetry enhances electron interactions near 
half-filling as one would expect from the increase in the 
density of states (46). 
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In what follows we consider the case of a very strong 
magnetic field <p > 1. The latter condition corresponds 
to an exponentially large effect of the field on the single- 
electron NT spectrum. For that reason we will neglect 
the Zeeman effect, which is linear in field. 

We consider the interacting problem whose Hamilto- 
nian in the forward scattering approximation 27 reads 



Htot — T~to + Hint 



(76) 



where 



and 



(77) 
(78) 



Here the Hamiltonian (77) describes four noninteracting 
fcrmion flavors (4 = 2 sp ; n x 2 va iiey) in a nanotube with the 
bare gap A = hvS/R, subject to a perpendicular mag- 
netic field, ip(y) — 2<j)sm{y/R), ip' y = d(p/dy, where the 
dimcnsionlcss field strength <\> is defined in Eq. (16). The 
Dirac spinors ^ a are operators in the second-quantized 
representation. The Coulomb interaction between elec- 
trons is described by the Hamiltonian (78), where the to- 
tal density in the forward scattering approximation reads 



circumferential degree of freedom using the following sep- 
aration of scales. The effects of magnetic field occur on 
the short scale of the tube radius R, for which the rele- 
vant energy scale is ~ fiv / R. The effect of the Coulomb 
interaction between electrons accumulates over a length 
scale that is much greater than R, as described below. 
Therefore the effective ID description of the interacting 
NT electrons can be obtained by first integrating out the 
circumferential coordinate y in Eq. (76) and then taking 
into account the Coulomb effects. Thus we obtain the 
effective one-dimensional Hamiltonian 



J a 



(83) 



with the bare velocity reduced due to supersymmetry, 



5o (</>) 

similar to Eq. (44). Here 50 is given by Eq. (51), 



_ ( Xlc 
X2c 



(84) 



(85) 



the one-dimensional electron density (calculated from 
half-filling) 



(79) 



XaXo 



(86) 



with the 2K n harmonics [K defined in Eq. (3)] neglected. 
The electron-electron interaction potential in the pres- 
ence of a substrate with a dielectric constant e is 



V(x,y) = ^V (x,y) , 



where 



Vo(x,y) = 



Vh 2 + x 2 



h = 2Rsm{y/2R) 



(80) 



(81) 



The problem (76) is SU(4) invariant with respect to ro- 
tations in the space of the four fermion flavors 

Below we focus on the low-energy properties of the 
problem (76). This allows us to utilize the projection on 
the supersymmetric basis (39): 



*«W=»a {x)^f ] (y) + X 2 a (x)^> (y) 



(0), 



(82) 



where the factorization of motion along x and y holds 
due to the assumption that the magnetic field does not 
vary along the tube. Using Eq. (82) we will reduce the 
problem (76) to the one-dimensional one, bosonize it, and 
estimate the effect of interactions on the plasmon velocity 
and on the semimetallic gap. 

Let us perform a projection of the problem (76) onto 
the basis (82). This can be done by integrating out the 



and the ID interaction potential at kR 1 



2e 



V(k)^— -ln[l + (fci?)- 2 ] 

£ + 1 



(87) 



In writing Eq. (83) we dropped the terms of the order 
[e 2 /(e + l)] In that are small compared to V(k) at kR <C 
1. These terms appear since at <fi > 1 the states \ia 
and X2 a are localized on the opposite sides of the tube. 
Therefore strictly speaking the interaction between the 
same components of the spinor (85) is cut off on the scale 
of magnetic length (64) rather than of the NT radius. 

With the difference between the short-distance cutoffs 
in the potential (87) neglected, the effective ID Hamilto- 
nian (83) remains SU(4) invariant. It can be bosonized 



in the standard way, Xa oc e a . This procedure im- 
mediately yields the renormalized plasmon velocity for a 
metallic nanotube, 



e = *va ( 0_L- 
90(9) 



(88) 



where the charge stiffness (or the dimensionless interac- 
tion strength) 



KM = 1 + 



TXhV 
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is enhanced by the magnetic field. (By the tilde we de- 
note the physical quantities in the presence of ID inter- 
actions.) Thus the plasmon velocity suppression factor 
go(<fi) = [v(<f>) /-5(0)] _1 due to the magnetic field is given 

by 



9o(<f>) = 



K(0) 



1 1/2 



K{4>) 



(90) 



It is reduced compared to the noninteracting value go (</>) 
because of the enhancement of the interaction strength 
due to the perpendicular magnetic field. For a large inter- 
action strength K 1, the effect of electron interactions 
on the supersymmetry is dramatic: 



9o(4>) - [go(<j>)} 



1/2 oc 



(91) 



effectively reducing the field strength <j> > 1 by a factor 
of 2 in the exponential. 

Consider now the semimetallic gap in the presence of 
a magnetic field. In bosonized language, this gap is esti- 
mated as the energy of a composite soliton of the charge 
and flavor modes, 26,28 with its energy dominated by that 
of the charge mode at K ^> 1. The essential feature 
for the present analysis is that the effect of the mag- 
netic field on the Gaussian part of H e s factorizes, renor- 
malizing the velocity v, with the backscattering term 
8/R = A(c/))go((j))/fiv inside the integral in Eq. (83) inde- 
pendent of the held. A straightforward calculation shows 
that the perpendicular field reduces the renormalized gap 
by the factor (90) obtained for the plasmon velocity: 



A(0) 



9o{4>) 



A(0) ~ K 1/2 (0)L> 1/5 A 4/5 (0) . (92) 



Here A(0) is the (renormalized) semimetallic gap in the 
absence of the field and D ~ TivjR is the one-dimensional 
bandwidth. 

In Eq. (92) the value of K is assumed to be taken at 
ql c h ~ 1, where Z c h 3> R is the size of the charged soliton 
in the bosonized description. 26 The universal power law 
4/5 in the gap renormalization (92) is valid in the limit 
K 1. Equations (88) and (92) show that the character- 
istic supersymmetry features of the velocity and minigap 



suppression persist in the presence of Coulomb interac- 
tions. However, the effect of the perpendicular magnetic 
field is strongly reduced by the electron interactions due 
to the density-of-states increase (46). 



VIII. CONCLUSIONS 

To conclude, we have shown that the interesting prop- 
erties of the nanotube electron spectrum in a perpendic- 
ular magnetic field found in Ref. 7 can be understood as 
a consequence of supersymmetry of the low-energy NT 
Hamiltonian. We have demonstrated that supersymme- 
try ensures stability of the zero-energy state in metallic 
NT's and yields a corresponding enhancement of the den- 
sity of states. In semiconducting NT's, supersymmetry 
leads to an energy gap suppression that can be observed 
in transport or tunneling measurements. We also consid- 
ered the effects due to the trigonal warping interaction 
arising from higher-order gradient expansion terms that 
violate supersymmetry and lead to field-sensitive mini- 
gaps in the metallic NT spectrum. Finally, we have found 
that supersymmetry persists in the presence of electron 
interactions, but the reduction of both the renormalized 
plasmon velocity and the excitation gap is weakened due 
to effectively increased interaction strength. 
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